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We show that the formal asymptotic limit for all rotating polytropic axisym- 
metric perfect MHD flows is a kinetic energy dominated wind which collimates 
to paraboloids around the symmetry axis. We reach this result by showing that 
another, a priori possible, solution with finite Poynting flux can be excluded on 
the following physical grounds: (1) the wind velocity does not exceed the fast 
mode speed everywhere and (2) the circumpolar current increases with distance 
5h ' from the source. 

We show that asymptotic hoop stress collimation is mathematically robust 
and we give strong arguments why recent 'anti-collimation' claims are not correct. 

However, in practice, due to the very slow logarithmic decline of the circum- 
polar current with increasing distance from the source, there is a broad inter- 
mediate regime with significant Poynting flux. This intermediate regime, rather 
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than the mathematically exact asymptotic regime, may well apply along the fi- 
nite length of the jet. We briefly discuss peculiarities that would be associated 
with Poynting jets in the intermediate regime. Force-free initial conditions in the 
near field are most likely to produce such jets, in which most of the energy flux 
is electromagnetic. 



1. Introduction 

We have previously established (Heyvaerts and Norman (1989)) that any stationary 
polytropic axisymmetric magnetized wind will collimate to paraboloids or cylinders along the 
symmetry axis at large distances from the source according to whether the electric current 
brought to infinity per hemisphere by the wind is finite or vanishes. In two companion 
papers (Heyvaerts and Norman 2002a,b) we have presented explicit asymptotically-matched 
solutions for both classical and relativistic winds. Our results however have left open the 
question of whether the asymptotic circumpolar current vanishes or not. We now address 
this issue which we already touched in earlier reports on our work (Heyvaerts and Norman 
(1997), Heyvaerts (1999)). 

Previous work has described the transfield equation (Okamoto (1975) Heinemann and 
Olbert (1978)) that expresses the balance of forces in a general magnetized rotating flow in 
the presence of a gravitational field. Blandford and Payne (1982) showed how collimated 
flows from disks could be obtained with a similarity solution associated with a particular scale 
free current and field distribution on a Keplerian disk. Contopoulos and Lovelace (1994) have 
found a similarity solution incorporating gravity and pressure, extending and generalizing 
the Blandford and Payne approach, as also did Ostriker (1997). Lovelace and collaborators 
have also considered relativistic jets, using the force free approximation (Lovelace (1976), 
Lovelace et al. (1991), Lovelace et al. (1993)). In Lovelace et al. (2002) they find collimated 
Poynting flux jets with surrounding winds. Numerical solutions of the steady state rotating 
axisymmetric MHD system have been presented by Sakurai (1985) and Sakurai (1987) for 
the split monopole and magnetized disk. A recent study of the relativistic case is given by 
Bogovalov (2001). Collimation along the axis is clearly evident in both cases. This result 
has been strengthened by a number of analytical and numerical studies (Sauty et al. (1992a), 
Sauty et al. (1992b), Sauty and Tsinganos (1994)). Shu and collaborators, have analysed 
protostellar outflows (Shu et al. (1994a), Shu et al. (1994b), Najita and Shu (1994), 
Shu et al. (1995)). Pelletier and Pudritz (1992) produced special solutions that exhibited 
focussing and recollimation. Numerical results (Ustyugova et al. (1995), Ustyugova et al. 
(2000), Ouyed and Pudritz (1997a), Ouyed and Pudritz (1997b), Krasnopolsky et al. (1999)) 
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further confirmed this. 

Heyvaerts and Norman (1989) have used asymptotic analysis to understand the general 
properties of the shape of the field lines and the flow structure at large distances from the 
source as a function of the conserved flow quantities. They have shown that the structure of 
the field and flow at infinity is controlled by the total poloidal electric current flowing in the 
wind about the polar axis. In the case where this current vanishes, the poloidal field (and 
flow) lines asymptotically approach to parabolae which focus to the polar axis in such a way 
that both r — > oo and (z/r) — > oo where standard cylindrical coordinates (r,9,z) are used 
with z along the symmetry axis. An expression giving the shape of these lines at infinity 
far from the polar axis, as a function of the flux function was derived. When the circum- 
polar poloidal current is non-zero, the poloidal field (and flow) surfaces have been shown 
(Heyvaerts and Norman (1989), Heyvaerts and Norman (1996)) to asymptotically approach 
to cylinders (possibly nested in asymptotically conical flux surfaces). The asymptotic radii 
of these cylindrical flux surfaces are given as a function of magnetic flux in the super- Alfvenic 
regions. In a relativistic generalization of our 1989 paper, Chiueh, Li and Begelman (1990) 
showed our results held not only for the non-relativistic case but also for the relativistic case. 

In a recent paper Okamoto (2003) insists that the asymptotic circum-polar current must 
vanish and that the asymptotic structure of magnetic surfaces must be of the conical type. 
Our earlier results (Heyvaerts and Norman 1989) by no means imply that the asymptotic 
circum-polar current should assume a non- vanishing value, contrary to what Okamoto (2003) 
claims. This point has been made clearly in a number of conference reports (Heyvaerts and 
Norman (1997), Heyvaerts (1999)). On the contrary, our present results, which extend 
and make precise previous arguments published in these reports, just conclude the contrary, 
i.e., that the circumpolar current vanishes aymptotically in the mathematical sense. Thus, 
there is no disagreement between us on this point. By contrast, Okamoto's antifocusing 
statement is inconsistent with an asymptotically vanishing circumpolar current. A non 
focusing magnetic surface must obviously be a surface such that (z/r) approaches a finite 
value at infinite distance from the origin on that surface, or, if such a limit does not exist, 
it is a surface along which (z/r) remains bounded from above. On such a magnetic surface 
pr 2 must approach a finite value (or remain bounded from below) at infinite distance from 
the origin. Indeed, an elementary calculation shows that when (z/r) approaches a finite 
limit (or is bounded) so does r|Va|, which implies that the polodal current enclosed in 
this surface approaches a finite value (see Heyvaerts and Norman (1989) and Heyvaerts and 
Norman (1996) for details). If the mathematical asymptotic state, because of a very slow 
convergence, is to be reached only at distances larger than the actual length of the jet, there 
is of course no point to confront the actual wind structure with this asymptotic solution. 
While this is the point made in this paper, it is not what is meant by Okamoto (2003) when 
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he claims conical asymptotics. 

The aim of this paper is precisely to determine whether, in general, the asymptotic 
circumpolar current vanishes or is finite. This question cannot be dealt with in the framework 
of a self-similar model, nor of any model which imposes a priori constraints on the solution. 
In this paper we discuss this issue in terms of the properties of the first integrals of the 
motion, which we take as given. These first integrals are determined by boundary conditions 
and by the criticality conditions. We consider both classical and relativistic winds. 

The paper is structured as follows. In Section 2 we review the basics of both relativistic 
and non-relativistic stationnary axisymmetric rotating MHD winds. In section 3 we summa- 
rize results on the asymptotic structure of MHD winds relevant to the distribution of flux 
in the asymptotic domain. In section 4 we show that a necessary condition for solutions 
with a non-vanishing circumpolar asymptotic current to be possible is that the function 
a(a)E(a)/Q(a) has a minimum value, I sup , away from the polar axis. We show that this 
non-vanishing asymptotic circumpolar poloidal electric current is proportional to this mini- 
mum value. The analysis extends to relativistic winds, for which a non-vanishing asymptotic 
proper current must be the minimum value, K sup , of a(a)(E(a) — c 2 )/f2(a). We next show 
in section 5 that classical winds which are super fast mode at infinity must have a current 
less than ^I sup - A similar upper bound is derived for relativistic winds. Due to the very slow 
decline of the current, the final asymptotic regime may in practice not be reached, allowing 
for intermediate quasi-asymptotic regimes. These are discussed in section 6 where we em- 
phasize systems carrying close to the maximum current. These winds split into a separate 
conical jet of small opening angle and a circum-equatorial wind separated by a region where 
the energy flux is almost all in Poynting form. We summarize our results in section 7. 

2. General Properties of MHD Winds 
2.1. Notation and Definitions 

We study perfect MHD, axisymmetric, stationary, polytropic jets and winds. The physi- 

— * 

cal quantities are denoted by their usual symbols, p, p, v, B for the density, pressure, velocity 
and magnetic field respectively. The gravitational potential is $g- Any axisymmetric vector 
quantity can be split into poloidal and toroidal components. The magnetic flux $ m through 
a circle of radius r centered on the axis at altitude z is written as $ m = 2na(r,z). The 
poloidal part of the magnetic field can then be expressed as 

-» 1 da _ 1 da _ 

Bp = — —e r + - —e z (1) 

r oz r or 
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The magnetic surfaces, generated by the rotation of field lines about the axis, are surfaces 
of constant a(r, z). The value of a suitably labels them. 



2.2. First Integrals 

A rotating, axisymmetric, stationary, polytropic perfect MHD wind flow admits five 
first integrals of the motion (see for example Heyvaerts (1996)) a, f2, L, E and Q which are 
conserved on magnetic surfaces, a(r,z). For classical winds, these first integrals are defined 
by: 

pv P = a(a)B P , (2) 
pv e = a(a)B e + prVL(a) (3) 

L ^ = rv '-£ty (4 > 

2 2 1 — 1 p p a{a) 

V = Q(a)p r (6) 
T is the polytropic index. The Alfven radius r^(a) and the Alfven density pa(o) are defined 

by 

r 2 A (a) = L(a)/Q(a) p A (a) = p a 2 (a) (7) 

The Alfven point on a magnetic surface a is at r = r^(a). The density at this point is 
Pa- The toroidal variables Vg and B 9 can be expressed in terms of first integrals and of the 
density p as: 

L p L-r 2 Q 

ve = - + - 5 , (8) 

r r p a z — p 

p ap L-r 2 Q 

B e = 2 • ^ 

r p Q a z - p 



2.3. Relativistic Winds 

Relativistic flows are characterized by a local Lorentz factor, 7, defined by: 

do) 
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The proper rest mass density is denoted by p and the specific momentum (the momentum 
per unit mass) is 

U — ^fV (11) 

The proper gas pressure is assumed to be related to the proper density by Eq.(6), where Q 
is constant following the fluid motion. We define the function 



which is also equal to (1 + J dP/pc 2 ) calculated at constant entropy for a polytropic gas. 
We denote the electric charge density by p e , the electric current density by j and the mass 
of the central object by M*. The special-relativistic equation of motion can be written as: 

- / CM. \ 

1P {v • V)( 7 ^) = "VP +jxB + Pe E + 7P V f 7 £-^ J (13) 

The relativistic form of the laws of mass conservation, isorotation, angular momentum con- 
servation and Bernoulli involve surface functions E, a, L, Q and Q, which are in this case 
defined by the relations: 

ipvp = a(a)Bp, (14) 
1 (v e ~rn(a))=a(a)B e /p, (15) 

7 ^-^) =L(0) ' (16) 

2 GM rn(a)B 

7Uc 5-) ^ = E ( a ) 17 ) 

R poa(a) 

Note that E in Eq.(17), now includes the rest mass energy. The rotation rate of the magnetic 
field, Q(a), which appears in equations (15) and (17), is defined in terms of the electric field 
by: 

E = -n(a)Va (18) 
Toroidal variables may be expressed using Eqs.(16) and (15): 

rB e =p ap (19) 



7^, = ? + Z ^T^F^ (20) 



L + p L_ Z ir 2 & 
r r PqQ 2 ^ — p 

Since 7 depends on vg, the elimination of the toroidal variables in Eqs.(19) and (20) is not 
yet complete. These expressions can be substituted in the Bernoulli equation (17), which 
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can then be solved to obtain an expression of 7^ in terms of the poloidal variables. This 
eventually gives the toroidal variables in terms of poloidal variables alone as: 

_ L^-GMJR)-t^E 

rBe ~ _//oC V - GM,/i2)(l - tooKIp) - (21) 

n* = *l(l ^-g^/g-^ ) (22) 
\ p (1 - /io« ^/ p)r 2 ilE - r 2 ilL ) 

These expressions have to be regular where their denominator vanishes, defining quantities 

which play the role of Alfven density and Alfven radius in the relativistic context. 



3. Asymptotics of MHD Winds 
3.1. Semantics 

We define important expressions to be used in this paper. The asymptotic domain 
consist of all points which are located on their own magnetic surface far away from the 
corresponding Alfven point. A magnetic surface is said to be asymptotically parabolic if the 
limits following this surface of r and (z/r) are both infinite. A magnetic surface is said to be 
asymptotically conical if r approches infinity and (z/r) approaches a finite limit, tan^ 00 (a). 
This does not imply that (z — rtanip^) approaches a finite limit: conical magnetic surfaces 
may have parabolic branches. An asymptotically cylindrical magnetic surface is one on which 
r approaches a finite value, r 00(a). The ratio r/r^ then also approaches a finite limit. When 
r becomes infinite on a magnetic surface, the surface is said to flare out. 

A neutral or null magnetic surface is one on which the poloidal field vanishes. The 
toroidal field also vanishes, since it is generated from the poloidal field by rotation. A neutral 
surface is strictly speaking not a magnetic surface, since it contains no magnetic field lines. 
Null surfaces are sandwiched between regular magnetic surfaces, so that the definition still 
retains meaning in the limit. The immediate vicinity of neutral magnetic surfaces is of special 
interest since they are exceptional regions where electric currents can flow in the asymptotic 
domain (Heyvaerts and Norman (2002a), Heyvaerts and Norman (2002b)). We refer to these 
regions as neutral surface boundary layers. The polar boundary layer in the vicinity of the 
polar axis is similarly a region where electric currents can flow in the asymptotic domain. 
We call the region away from boundary layers the field region. 

The total poloidal current, J = 2itrB e / 'p 0: is the electric current through a circle of axis 
z passing the point (r, z). For positive fl and a, J is negative. We shall loosely refer to the 
quantity 

/ = -rBe/p (23) 
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as the current. The total poloidal electric current enclosed in a neutral surface is zero. This 
means that neutral surfaces separate the wind in a number of cells in each of which the total 
current separately closes. We refer to the space between two neighbouring neutral magnetic 
surfaces as being a cell. In the asymptotic domain, current-carrying regions are restricted to 
regions of small extent (Heyvaerts and Norman (2002a), Hey vaerts and Norman (2002b)). 

There is a total current that flows about the polar axis which we refer to as the circum- 
polar current. This circumpolar current varies with the distance to the wind source. 



3.2. Relevant Results 

For r ^> tai the azimuthal velocity vanishes while the current approaches the value 

pr 2 Q(a) 



Hoa(a) 



(24) 



Heyvaerts and Norman (1989) show that pr 2 is bounded from above. The toroidal component 
of the velocity, vg, approaches zero as r approaches infinity on any flaring magnetic surface. 
Thus the flow velocity becomes poloidal while / remains bounded. If (pr 2 ) approaches a finite 
limit, the asymptotic structure of the magnetic surfaces consists of a set of conical surfaces, 
inside of which cylindrical magnetic surfaces are nested. Such flows convey Poynting flux to 
infinity. On cylindrical magnetic surfaces none of these results strictly apply. Here, as in 
our previous papers, we nevertheless assume that the asymptotic cylindrical radius on any 
magnetic surface is much larger than the corresponding Alfven radius. If pr 2 does converge 
to zero at large distances, the asymptotic structure of magnetic surfaces consists of nested 
paraboloids (Heyvaerts and Norman 1989). No Poynting flux reaches infinity and all the 
energy emerges in kinetic form. 

A detailed asymptotic solution in terms of given first-integrals has been presented in 
the companion papers (Heyvaerts and Norman 2002a,b) for both classical and relativistic 
winds. We have shown that electric currents are distributed in the asymptotic domain in 
boundary layers about the pole and in the vicinity of neutral magnetic surfaces, with very 
little current density flowing outside these regions. Detailed solutions for the inner structure 
of these current-carrying boundary layers were also discussed. Let a n be the flux variable of 
some neutral magnetic surface. If there is a mass flux along the surface, the first integral 
a diverges for a — > a n , as can be seen from Eq.(2). We have established that outside these 
boundary layers the transfield equation takes the simple form 

/ = /«>(&) (25) 
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where b labels surfaces orthogonal to magnetic surfaces. Specifically, the label b is taken to 
be the value of the ^-coordinate on this orthogonal trajectory at the polar axis. Variables 
a and b could, in principle, be taken as space coordinates replacing r and z, with a playing 
the role of an angular variable and b the role of a radial variable. 

The current variable, J, in general depends on both a and b. Eq.(25) expresses the fact 
that its dependance on a almost completely disappears in the field-regions of the asymp- 
totic domain, outside the current carrying boundary layers. More precisely, /<x>(&) is only 
piecewise constant, reversing sign when crossing from one cell to another. This follows from 
consideration of the equilibrium of the neutral sheet boundary layers where the dominant 
toroidal magnetic pressure on each side must balance the other. So, Bg, and / reverse sign 
as neutral magnetic surfaces are crossed. This means that the absolute value of the total 
poloidal current |ioo(&)| passes unchanged from one to the next cell. 

Winds therefore could come in two separate classes: Poynting jets approaching a finite 
value of \Ioo(b)\ or kinetic winds having asymptotically vanishing current. Whether both 
classes actually occur in nature remains to be discussed. 

Outside the current-carrying boundary layers, the Lorentz force j x B vanishes since 
jp and Bp become both small. Inside boundary layers, other forces balance the Lorentz 
force (Sauty et al. 1999). The gas pressure is the most likely additional force in the case of 
polytropic flows since its decline with r is the slowest. The polar boundary layer has the 
structure of a column pinch, while the neutral surface boundary layers have the structure of 
sheet pinches (Heyvaerts and Norman (2002a)). 

3.3. Distribution of Flux on Orthogonal Trajectories: Classical Winds 

For large values of r/r^, it is possible to determine the distribution of flux along tra- 
jectories orthogonal to magnetic surfaces by quadrature. Denoting the terminal velocity on 
surface a by i>oo(a), equations (2), (5) and (9) simplify to the form: 



Eliminating and p we obtain an expression for r| Va| in terms of the first integrals and /: 




(26) 




(27) 



(28) 




(29) 
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Let the curvilinear abcissa along the trajectories orthogonal to the magnetic surfaces be 
denoted by a, conventionally increasing from pole to equator. Then |Va| = da/da, and the 
relation (29) becomes: 

da_ = n(a)da 

r Haifa, b)V2y/E - /(a, 6)fi(a)/a(a) 

This relation holds true both in the field-regions, where I(a, b) does not depend on a, and 
in current carrying boundary layers, where this dependence is rather strong. 

The actual form of the quadrature relation (30) depends on the actual shape of the 
orthogonal trajectories to magnetic surfaces. These trajectories have been shown to be ap- 
proximately circles when I(a, b) approaches a non- vanishing limit 1^ (Heyvaerts and Norman 
(2002a), Heyvaerts and Norman (2002b)). This result has been extended in a WKB sense to 
the case when I slowly varies with b. Asymptotically cylindrical surfaces are approximately 
orthogonal to circles too. The position on such a circular orthogonal trajectory is specified 
by a latitude angle ip from the equator, so that 

z = r tan(^(a)) (31) 

For approximately circular orthogonal trajectories Eq.(30) integrates to 

tan ^ o = smh / - } == + tan ^ fll (32) 

\Ja V2noI{a, b) y/E(a') - I(a, b)tt(a')/a(a') J 

where a\ is some reference flux. Along the orthogonal trajectory b the current I(a, b) vanishes 
near the polar axis proportionally to a and causes tan(^(a, b)) to diverge as a approaches 
zero. Outside boundary layers, I(a,b) becomes a function Ioo(b), essentially independent of 
a at constant b, and the solution (32) simplifies. It will be shown below that the integral in 
Eq.(32) converges at a neutral surface. If, moreover, we neglect the flux in neutral boundary 
layers the integration in Eq.(32) can be extended from the equator, at flux a = A, where 
tanip(A) = 0, to the surface a with a constant value of \I(a,b)\ = 1^. This gives an 
approximate expression for ip(a): 

tan (('(»)) = sink I / — ^ ^(a')da' \ ^ 
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3.4. Distribution of Flux on Orthogonal Trajectories: Relativistic Winds. 

A similar analysis can be done in the case of relativistic winds (Heyvaerts and Norman 
(2002b)). In the field- regions the quantity K, defined by: 



K = ^-y/c 2 B 2 -Q 2 \Va\ 2 (34) 

becomes a constant. We refer to K as the proper total current (Heyvaerts and Norman 
2002b). K is r times the asymptotic value of the electromagnetic invariant (B 2 — E 2 /c 2 ) 1 / 2 . 
It reduces in the classical limit to the poloidal current / (see Eq.(50)). In field-regions K 
becomes independent of a along any orthogonal trajectory, giving: 

K = K(b) (35) 

The asymptotic relativistic Bernoulli equation can be written as: 

V^ 2 ~ (c 2 + KQ/a) 2 
u Poo = Too.Poo = c (c2 - KQ/a) (36) 

Note that in relativistic dynamics E contains the rest-mass specific energy, c 2 , so that E c i ass = 
Ereiat — c 2 ■ The relativistic form of Eq.(30) is: 

da = Q(a)^c 2 -t&(a) da 
r n cK(a,b) v^a) 

For relativistic winds, as for classical ones, the orthogonal trajectories to magnetic surfaces 
at large distances closely resemble circles, in which case Eq.(37) can be integrated to find 
the latitude -0(a) of the magnetic surface a on the orthogonal trajectory b. This gives: 



. , / f ai 1 \!(a- )<ia-{<- i l\{a.l))\!(a- )/<\{a' )) \ 

tan(^(a,6)) = sinh / ^ ^ 1 y y ' ' y " y " \ +tamp(a u b) 

\Ja ^ocK(a,> 



tt(a')da'(c 2 + K{a, b)tt(a') / 'a(a')) 

6) \]E 2 (a') - (c 2 + K(a,b)n(a')/a(a')) 2 

(38) 

where a\ is a reference flux in the same cell as a itself. In the case of dipolar symmetry, only 
one cell would be present per hemisphere. Neglecting the flux in the equatorial boundary 
layer, a\ is the equatorial flux, A, and tan^A, b) vanishes. Outside boundary layers, K(a, b) 
becomes a constant K^b), independent of a and the solution given by Eq.(38) simplifies as: 

, . ." A 1 n(a')da'(c 2 + K 00 (b)Q(a')/a(a')) . 

tan (#»,&)) = sinh | / - ±J K °° U UM = | (39) 

,a ^ oC Kco[D) ^E 2 {a') - (c 2 + K 00 (6)fi(a')/a(a'))" 
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4. The Circumpolar Current 

4.1. Mixed Cylindrical-Conical Asymptotics 

When I (a, b) converges to zero the solution, which consists of nested paraboloids, fills 
all space. We have shown in Heyvaerts and Norman (2002a) that, conversely, if |J(a, b)\ 
approaches a non-vanishing value at infinity, the asymptotic structure of magnetic surfaces 
consists of cylindrical surfaces, possibly nested in conical ones. These two geometries seem 
to clash because their common existence implies the presence of a big spatial region with 
no magnetic surface, between the last cylindrical surface and the first conical surfaces. A 
smooth connection is nevertheless possible. Let a* be the flux at which the magnetic surfaces 
switch from being cylindrical to being conical. The solution is continuous if r^a) approaches 
an infinite limit as a — > a* and if ip(a) approaches n/2, as a approaches a* from above. Thus, 
the difficulty of matching cylindrical to conical asymptotics is not sufficient to eliminate the 
idea that the total circumpolar current could remain finite. The relevant question is rather 
whether the wind's physics allows this. In a cylindrical region, r^a) cannot approach infinity 
uniformly as a approaches a*. Indeed, the limiting magnetic surface has the character of 
a paraboloidal surface described by z = F*(r) . Any asymptotically cylindrical magnetic 
surface with terminal radius r^a) is nested inside z = F^(r 00 (a)). Since F*(r) diverges with 
r , the limit r^a) can obviously not be reached uniformly. For future reference let us note 
that the terminal radius, r^a) is given by: 

roc a )=nexp / -= V ; 40 



4.2. Asymptotic Current for Classical Winds 

Eqs.(32) and (38) give the flux distribution explicitly in the field regions. Separate 
solutions have been obtained in current-carrying boundary layers and matched to the field- 
region solutions. For classical winds, a smooth asymptotic matching implies a relation, which 
we refered to as the Bennet relation, between the current 1(b) and the axial density at the 
distance b, po(b): 

p — 7 QoPo (b) = (41) 

1—1 Cto 

Matching also imposes the relation: 

= ( 2 -r)>n(M*„ + ln(f) m 
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Here no (b) = po(b)/ Pao is a dimensionless measure of the axial density and p^o is the Alfven 
density on the polar field line. The reference length £ is defined by 



T QoPao 1 

and A is the following integral, which depends on 1(b), that is, by Eq.(41), on n (b) or po(b) 



1 - T^VW- < 43 > 



Jo n s/E{a') - I(b)Q(a')/a(a') a 

As b approaches infinity, some term of Eq.(42) has to balance the divergence of the log- 
arithmic term on the right. If 1(b) asymptotically vanishes, so does po(^)(Eq.(41)), and 
A approaches a finite value. In this case, Eq.(42) is satisfied with po(b) approaching zero, 
consistent with Eq.(41). If, by contrast, 1(b) approaches a finite value, 1^, so does n (b) 
(Eq.(41)). The divergence of the second logarithmic term on the right hand side of Eq.(42) 
can then only be matched by a divergence of \(n (b)). As discussed in Heyvaerts and Nor- 
man (2002a), this implies that, as b approaches infinity, the circumpolar poloidal current 
1(b) rises to its maximum allowed value. This value is the absolute minimum of the function 
(aE/fl), reached away from the polar axis at a* (see Fig.l). Then: 

it- t ^ r . ,a(a)E(a). . , . 

I /oo I = hup = Mm^o I y (Q | ; I (45) 

For this current both tan(^ 00 (a)) and the terminal radius r^a) would simultaneously diverge 
as a* is approached. 

From Eq.(40) we see that \I\ should in this case increase with b. However, this result 
gives a direct contradiction since the circumpolar current must close back to the source 
via a low current density flow in the field-region. This picture implies that \I(b)\ must 
approach its limit from above. A decreasing \I(b) \ accelerates the outflow by the coiled spring 
force associated with the gradient of the toroidal magnetic field pressure. A monotonically 
decreasing \I(b)\ corresponds to continued, though weakening, acceleration from the wind 
source to infinity. An increasing \I(b)\ is associated with a wind flow undergoing continuous 
deceleration and thus the Poynting flux increases with distance. Therefore we conclude on 
these physical grounds that the circumpolar current cannot remain finite at infinity. 
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4.3. Asymptotic Proper Current for Relativistic Winds. 

For relativistic winds, an entirely similar analysis (Heyvaerts and Norman 2002b) has 
shown that relations similar to Eq.(41) and (42) hold true, namely 

QoPo (o) = (46) 



r — i o 
^ = P -r, M „„, + ,g'| ,4 7) 



where X r , defined by 



A = f A da' Sl(a!)y/<? + K(b)n(a')/ a (a') V E£ - (c 2 + K(b)Q /a o y 
r Jo «o y/E 2 (a') - (c 2 + K(b)n(a')/a(a')) 2 {l ^/c 2 + K(b)n /a 

is a function of n (b), since K(b) is given by Eq.(46). Similar arguments then show that, if 
K is to approach a finite limit, its value should be the absolute minimum of the function 
(a(E — c 2 )/Q), reached at a non-zero regime change flux, a*. So, if Koo is not zero: 

.a(a)(E(a) - c 2 ) , . . 

I^ool = K sup = Min a ^o | 1 -| (49) 

This limit, however, can only be approached from below, again leading to a physical contra- 
diction. Although K is not directly related to the poloidal current / anymore, it is possible 
to see that, if the wind source is the only current source, K should be a decreasing function 
of b. Indeed, the asymptotic limit of K, K^b), is related to the asymptotic current Ioo(a, b) 
enclosed in surface a at distance b by (Heyvaerts and Norman 2002b) 

/ooM) =loo(a,b)K oc (b) (50) 

where 7oo(a, b) is the Lorentz factor on surface a at distance b. Again, Joo(a,6) decreasing 
with b, at constant a, implies plasma acceleration, since the field-aligned poloidal Lorentz 

— * — * 

force j P x B e is then in the sense of the motion. This causes 7oo(a, b) to increase (Begelman 
and Li 1994). Thus, loo (a, b) decreasing with b implies that K^b) is also decreasing with b. 
Solutions with a finite asymptotic value of \K\, which can be approached only by increasing 
values of \K\, are therefore physically inconsistent. The circumpolar proper current cannot 
be finite at infinity for relativistic winds. 



4.4. Asymptotic Current: Progressive Deconfinement 



That the circumpolar current can only be zero or given by Eq.(45) can also be understood 
by the following analysis. Let us suppose that an external pressure independent of z confines 
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the jet in such a way that it has cylindrical geometry. At low confining pressure, its structure, 
which is entirely determined by the first integrals and by the boundary condition, will consist 
of a circumpolar current-carrying channel surrounded by a large, almost pressureless, field 
region carrying very little current density. In this latter region, the current enclosed in a 
magnetic surface a is almost independent of a and equal to /. This field-region value of / 
depends on the external confining pressure P ext . We calculate explicitly in Appendix (A) the 
relation between I and P ext , both when the edge of the jet is either a magnetic surface or a 
null surface. The conclusion of this analysis is that / approaches either zero or I sup (K sup in 
the relativistic case) as P ext decreases to zero. 

4.5. Asymptotic Current: Conclusion 

The discussions in subsections (4.2) and (4.3) have clearly shown that if the wind extends 
over very large distances, the amount of circumpolar electric current should approach zero 
as z approaches infinity. When this terminal regime is reached, all of the wind energy is in 
kinetic form. As shown in Heyvaerts and Norman (1989), Heyvaerts and Norman (2002a) 
and Heyvaerts and Norman (2002b), the magnetic surfaces are then paraboloidal. 

5. Characteristic MHD Speeds 

5.1. Terminal Speed vs Fast Mode Speed 

In the absence of external pressure, a wind should eventually achieve, on each field 
line, a terminal velocity exceeding all the MHD characteristic velocities. We now show 
that the fully asymptotic wind regime cannot satisfy this condition if the wind is to carry 
a finite circumpolar current. Specifically, super fast mode velocities cannot be reached in 
the vicinity of the regime-change surface a* in the field-region. As discussed above, the 
terminal velocity has to vanish at this surface and, for a classical wind, the circumpolar 
current should equal I sup , the minimum value of the function (aE/fl). Does the terminal 
velocity, though vanishing, remains larger than the fast mode speed? A necessary, though 
not sufficient, condition for the flow remaining super-fast-mode is that the fast mode velocity 
itself approaches zero as a approaches a*. 

The position and velocity at the fast critical point are obtained from the Bernoulli 
equation (5). Let us write it, for a given magnetic surface a, in the form £>(r, p) = 0. The 



-16- 



critical points and associated velocities are obtained by solving the system 

r d A = _I ± r d 3L _ r ^£ + fiV + ^ L ' ~ ^ V = o (51) 
dr 2 p 2 rfr r 2 (p^ — p) 2 

«9£ c?B\ v _ x pp\ (L - ^r 2 ) 2 

~ P 5p" = "^ + rQP + ~ (p A -p)3 =° (52) 
The velocity at the fast point is then, for classical winds : 

v) = TQp 1 *' 1 + ^ PP \ ( f 2 ~ r ^ (53) 
J r\p A - p) 6 



A lower bound to this is 



V > Pu 

f ~ r 2 (p A - p) 3 



(54) 



which, from Eqs. (9) and (2) means that vj > (Bp + £>|)/p p. A weaker, but accurate, 
lower bound for t> 2 is Bg/p p, which is obtained in the limit r ^> and p <C p^. Using 
Eq.(24), the inequality (54) reduces, in this limit, to 

v 2 f >\ — I 55) 

The fluid terminal velocity v^i^a), given by Eq.(28), cannot remain larger than IQ/a if 
is to approach zero as a approaches a*, since |/| is supposedly close to a non-zero value and 
a supposedly remains finite. Note however that this would not be so if a* were a neutral 
magnetic surface, where / vanishes and a diverges. A more general argument is as follows. 
Eq.(28) and the inequality (55) indicate that the wind can only be super fast mode if : 

In the field-regions of the asymptotic domain, |/| is a constant, 1^. \I\ should be less than 
I sup , the minimum value of \aE/Q\, otherwise, from Eq.(28), v 2 ^ would not be positive ev- 
erywhere. If it is required that the wind be everywhere super-fast mode at large distances, 
inequality (56) shows that the upper limit I sup cannot be reached. In particular, the cylindri- 
cal asymptotics regime, for which 1^ = I sup , is not consistent with the wind being everywhere 
super-fast-mode. Winds in an intermediate regime, which deliver a circumpolar current very 
close to the maximum I sup , should still be in the sub-fast mode regime. Unless the envi- 
ronment of the jet makes it possible that its velocity does not exceed the fast mode speed 
asymptotically, the electric current should eventually fall below ^I sup at very large distances. 
This argument again points to the conclusion that the completely asymptotic regime should 
have a vanishingly small circumpolar electric current. 
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5.2. Fast Point in the Relativistic Regime 

Similar conclusions are found in the case of relativistic winds. The location and velocity 
of the fast point can be found as in the classical case. Neglecting gravity and the gas entropy, 
the relativistic Bernoulli equation takes the form 



Assuming the shape of field lines is given, Eq.(14) turns Eq.(57) into an equation for u p , 
B(u p ,r) = 0. At critical points the differential of B vanishes. The position, 77, of the fast 
point is then given by : 

where s is the curvilinear abcissa along a field line. Eq.(58) can be satisfied at large distances 
in any particular geometry compatible with the Bernoulli equation (appendix (B)). For an 
initial monopolar field at the source of a highly magnetized, relativistic wind the fast point 
is at a large distance from the light cylinder (Beskin et al. 1998). The specific momentum 
at the fast point is obtained from the condition that dB/du p vanishes. The Lorentz factor 
associated with the poloidal velocity at the fast point is: 

7,/ = (f ) 3 (59) 

Eq.(57) can be solved when the flow is in the asymptotic domain, since pr 2 becomes a con- 
stant at large distances, resulting in Eq.(36) which relates the asymptotic specific momentum 
to the proper total current K. A condition for the flow to be super-fast-mode, similar to the 
inequality (56), can be derived for relativistic flows. The fast mode speed, Cf, for a diffuse 
highly magnetized medium where the Alfven speed ca = (B 2 / pop) 1 ^ 2 ) may exceed the speed 
of light, is given by 

111 

+ (60) 

The Alfven speed associated with the total field, ca, is larger than the the Alfven speed 
associated with the toroidal field, cao- This gives for c/ the following inequality: 

r 2 r 2 
c -t- c Ag 

The proper current K is related to / and 7 by Eq.(50) and can be expressed in terms of p 
by (Heyvaerts and Norman 2002b): 

pr 2 fl . . 

K = 62 

p a 
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Using these relations, the inequality (61) can be transformed into the following inequality 
for the specific momentum Uf = cfy(cf) associated with cf. 

u) > E 2 ^ 2 (63) 

A necessary condition for the relativistic flow to be super-fast-mode is that the terminal 
specific momentum of the wind flow, wpoo, given by Eq.(36), exceeds the lower bound to the 
asymptotic value of Uf, given by Eq.(63). This can be written as: 

Now, if the asymptotic flow regime is to have a finite value of the circumpolar proper current 
K, the latter must be the minimum value of a(E — c 2 )/Q, reached at some non-zero a*. At 
this particular value of a, the inequality (64) cannot be satisfied, since its left side vanishes 
whereas its right side is strictly positive. Expanding (64) for E close to c 2 gives again the 
inequality (56) for classical dynamics. Therefore, relativistic winds in an intermediate regime 
having a circumpolar current very close to the maximum I sup should still be in the sub-fast 
mode regime in some regions. Since the fast mode critical surface is unlikely to be found far 
in the asymptotic domain, this situation should be exceptional. By contrast, the completely 
asymptotic relativistic regime should have a vanishingly small circumpolar electric current. 



6. Intermediate Asymptotic Regime 

6.1. Is the Mathematical Asymptotic Regime ever Reached? 

In the fully asymptotic regime, the currents / and K should approach zero with increas- 
ing distance b to the wind source. However, our calculations have shown that they do so only 
very slowly, decreasing as the inverse of the logarithm of the distance to the source. This 
opens the question of whether, in fact, when the flow reaches the terminal shock separating 
it from the outer medium, / has indeed vanished. 

If not, the region in the wind cavity at r ^> ta would be in an intermediate asymptotic 
regime, i.e., one in which / is still finite, though smaller than I sup . For classical winds, the 
flux distribution would then be as given by Eq.(32), with the value of po(b) and 1(b) in the 
asymptotic domain given by Eqs.(41) and (42). Extreme focusing would be possible for I 
close to its maximum value, I sup , defined in Eq.(45). However such a wind would not yet 
have made its transition through the fast critical point. Similar considerations apply to 
relativistic winds. 
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Such an intermediate asymptotic state is associated with a finite value of I < I sup (or 
K < K sup ). It is appropriate to discuss the properties of winds which carry a current close 
to the maximum possible value, since these have strong focusing properties and may be 
directly related to ubiquitous jet phenomena. We show below that such winds have very 
special properties. For example, when the current is close to I sup , the wind splits into a 
focused jet and an equatorial wind with radial geometry. In between these two components, 
there is only a very weak flow, but a large Poynting flux. Such highly focused winds must 
be sub fast mode as we show below. 



6.2. Poynting-Flux Dominated Winds: Polar Jets and Equatorial Winds 

We now discuss the structure of rotating MHD winds which carry a circumpolar current 
close to the maximum possible value. This current is treated as remaining constant in the 
intermediate region as discussed above. The wind's structure is given by Eq.(32). Suppose 
that / is close to its maximum possible value, I sup , the minimum of the function (aE/fl), 
reached at some a* 7^ 0. It is assumed that (aE/Q) has such a minimum. In this case, 
the integral in Eq.(32) would almost diverge at a*. As a result, all the flux at a < a* is 
strongly focused about the axis. The focusing is still conical however because divergence is 
not exactly reached. The flux at a > a* is spread about the equator. The terminal flow 
velocity in the vicinity of a = a* almost vanishes, while the angular extent of this region is 
large. This structure is represented in Fig. (2). 

A wind close to the maximum current then appears to have a rather peculiar structure. 
It is split into two main flow regions, a circumpolar flow channeled in a cone of very small 
opening angle, and a conically spread equatorial flow. These two regions are separated by a 
region with very little flow, where the outflowing energy is almost entirely in Poynting form. 
This result has been also obtained in numerical simulations by Ustyugova et al. (2000) of 
the flow from a Keplerian disk threaded by a dipole-like magnetic field. 

Again, from the discussion of section (4), this regime cannot be the eventual asymptotic 
regime. The latter necessarily involves a vanishing circumpolar current. The solutions 
described above occur in the intermediate asymptotic regime , which may persist over large 
distances. That I or K are close to their maximum possible value, I sup or K sup , has however 
here the character of an extra and rather conjectural assumption. 

In the particular case when the function (aE/Q) reaches its minimum at a = 0, the 
situation regarding the focusing of the flow is different. The flux function a is proportional to 
r 2 near the polar axis. Symmetry about the axis implies that any physical quantity be even 
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in r, and reaches a flat extremum at r = 0. This however does not imply that functions of 
the flux a should behave similarly. On the contrary, a linear variation with a should be the 
rule. Therefore when a* = and the circumpolar current is supposedly close to a E /Q , 
there is no divergence of tan^a) (Eq.(32)), as a approaches zero, due to the vanishing of 
the square root in the denominator, because the latter usually has a simple zero. Most of 
the mass flux is then spread about the equator in a wind shaped like a conical fan. This 
structure is represented in Fig. (3). There is little flow about the axis, since the velocity, 
equal to y/2(E — IQ,/ a), is small in this region. Eq.(32) shows that the wind nevertheless 
fills all space. Indeed, ip(a) approaches n/2 as a approaches zero because / vanishes on the 
axis. The flow near the axis would in this case be weak, the energy being transported mainly 
in electromagnetic form in a circumpolar flux tube subtending negligible flux. 



6.3. Neutral Magnetic Surfaces do not Focus Winds 

In section (6.2) above, a wind carrying nearly maximal current in an intermediate asymp- 
totic regime has been considered. This situation induces a strong focusing of magnetic sur- 
faces. Note that at a neutral surface, where / also vanishes, Eq.(32) does not give rise to a 
divergence of tan?/>(a), since the integral in Eq.(32) is convergent (Heyvaerts and Norman 
2002a,b). Neutral surfaces do not force the magnetic surfaces nested in them to cylindrical 
shapes. It is shown in Appendix(A) that a neutral magnetic surface at the jet's edge does not 
cause the radius of a pressure-confined jet to diverge. We conclude that currents returning 
at neutral magnetic surfaces do not cause any general focusing of the structure. Okamoto 
(1999) shows that they rather have a defocusing effect. 

7. Conclusions on the Asymptotic Regime of MHD Winds 

The major conclusions of this paper are given below. 

1. If there were to be a finite circumpolar current at infinity then, for classical winds, 
this current should have a value equal to the minimum value of the function aE/fl. 
This minimum value, I sup , is the maximum allowed current. If this maximum current 
is achieved,then there exists a particular magnetic surface, a*, where aE/Q reaches its 
minimum and where the total energy is carried by Poynting flux alone. 

2. The limiting maximum current, I sup , can, however, only be reached from below. In 
other words, as the distance from the source increases, the circumpolar current would 
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have to increase. This cannot be achieved fro winds emanating from a finite, central 
source of magnetic flux. Thus, we conclude that, for non-pathological cases, the max- 
imal current cannot be reached from below. Therefore, the asymptotic circumpolar 
current must vanish. 

3. A winds carrying a cicumpolarcurrent, /, that is bounded from below such that, / > 
|/ snp , cannot have a terminal velocity everywhere in excess of the fast mode speed. 

4. The circumpolar current declines to zero only inversely logarithmically with distance. 
An important consequence of this very slow decline is that there exists an extensive 
intermediate asymptotic regime, where this current is still finite. Thus, significant 
Poynting flux can be carried over very large distances. The termination of the wind, 
for example at a terminal shock, may occur before this current has decayed to zero 
and the a parameter is still large 

Similar conclusions apply to relativistic winds: 

1. The total circum-polar proper current K is, for relativistic winds, less than K sup , the 
minimum value of the function a(E — c 2 )/Q. Any finite value of this proper current 
at infinity equals K sup . If this maximum current were to be achieved asymptotically, 
there would exist a particular magnetic surface, a*, where a(E — c 2 )/Q reaches its 
minimum value and where the total energy is carried by Poynting flux alone. 

2. The upper bound to the proper current, K sup , can only be reached from below. As 
in the non-relativistic case, it is not physically possible to reach this upper bound for 
winds emanating from a finite, central source of magnetic flux. Thus, the circum-polar 
proper current should formally approach zero at infinity. 

3. Winds that carry a proper current such that inequality (64) is not satisfied cannot 
have terminal velocities everywhere in excess of the fast mode speed. 

4. The circumpolar total proper current declines to zero only inversely logarithmically 
with distance. A consequence of this very slow decline is that there is an extensive 
intermediate asymptotic regime where this current is still finite. Significant Poynting 
flux can be carried over large distances, and the wind may terminate in a shock much 
before this current has decayed to zero. 

5. In the intermediate asymptotic regime, a classical or relativistic wind carrying a proper 
current close to the maximum allowed value I sup (or K sup ) consists of a conical wind 
of very small opening angle about the pole surrounded by a region in which most of 
the energy flux is in Poynting form. 
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6. Such close-to-maximum Poynting Flux carrying winds with polar jets and equatorial 
winds must be sub-fast mode. External boundary conditions may significantly influence 
their nature. 

7. For force-free initial conditions in the near field, the energy close to the source is all 
Poynting flux. Thus, relativistic jets emanating from a strongly magnetized base can 
maximize the output in Poynting flux 
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A. Unconfined Jets as limits of cylindrical pressure-confined jets 

We analyze the structure of a pressure-confined cylindrical jet with a given set of first- 
integrals as the confining pressure is reduced. For completeness, we consider both the case 
of a jet bordered by an ordinary magnetic surface and the case of a jet enclosed by a neutral 
magnetic surface. This is appropriate since the electric current organizes itself into closed 
cells and because the presence of a neutral magnetic surface at the jet's edge slightly changes 
its structure, favoring space filling. 

Where a uniform external pressure imposes cylindrical geometry for the magnetic sur- 
faces, the orthogonal trajectories are straight lines perpendicular to the rotation axis. Taking 
into account the smooth matching with the polar boundary layer the radius r(a) of the cylin- 
drical flux surface a is found to be: 

= j_ go^t; exp ( r vw* , (A . 1} 

T - 1 ^1 1 \J ^Iy/E(a') - m(a')/a( ■ ' 

The variables / and p are related by the Bennet relation 



a" 



, .QoP^ — (A-2) 
1—1 «o 

They both depend on the external confining pressure, P ex t, as does also the radius r* of 
the outermost magnetic surface. The boundary condition imposes equilibrium between the 
externally applied pressure and the total pressure at the edge of the jet. The latter is the 
sum of toroidal magnetic pressure and gas pressure since poloidal magnetic pressure can be 
neglected. If a* is simply an ordinary magnetic surface, gas pressure can also be neglected 
and the boundary condition becomes: 



ft- = f£ < A - 3 > 



where r* depends on / as given by Eq.(A-l). In this case, / and P ext are simply related by 
T ^ ^ faoQofa 1 ^ exp / r n(a>)da> \ = 



Itoo J \Jo V2fi W E ( a ') ~ Ito(a')/a{a') J Pext 

As P ex t decreases, the left hand side of Eq.(A-4) should diverge, which leaves at most two 
solutions for 1,1 = and, if (aE/Q) has a flat minimum, / = I sup (Heyvaerts and Norman 
(2002a), Heyvaerts and Norman (2002b), Lery et al. (1999)). The value I sup is reached by 
increasing I as P ext decreases. 
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Some complications arise when the surface a* is a neutral one. The jet is then bordered 
by a neutral surface. The proof that I = or I sup in this unconfined limit should be 
extended to this case. The outskirts of the pressure-confined jet then consist of half a sheet 
pinch boundary layer, in which the transfield equilibrium reduces approximately, n being 
the unit vector normal to magnetic surfaces, to: 

n ■ V (q P v + = (A-5) 



2p 

Be is given by Eq.(24). Integrating accross the boundary layer, we get, taking into account 
the divergence of the function a at the null surface: 

1 n 2 r 2 

Qp t + ^^ = Q*pI (A-6) 

2 p a 2 

where is the density at the jet's edge. It does not vanish, because of the confining pressure. 
Treating Q, E and Q as constants equal to Q*, -E* and f2* in this external boundary layer, 
Eq.(A-6) provides an expression for a as a function of the parameter 

X = — 77T (A-7) 
P*{b) 

namely: 

1 2Q^ 2 ( 1 



p a 2 Q 2 r 2 \X 2 X 2 - r J ^ A " 8 ^ 

Eq.(A-8) implicitly gives a as a function of X for known a(a). The function a cannot be 
treated as approximately constant because of its divergence at a = a*. Nevertheless, in our 
analysis, a is a given function of (a* — a) in this vicinity. It is expected that r itself should 
vary only little in the boundary layer if the latter is indeed thin, which will be checked 
a posteriori. The flux distribution versus radius results from Eq.(2), the velocity being 
obtained from Eq.(5). Disregarding the gravitational potential and the toroidal velocity and 
assuming that the velocity at the jet's edge is largely supersonic, we obtain 

dr a(a) 

r— = ^= A-9 

da X ^/2EZ 

Eq.(A-9) implicitly gives r(X) and Eq.(A-8) implicitly gives a(X). To derive an explicit 
parametric representation of the flux radius relation in the boundary layer we deduce a(X) 
from Eq.(A-8). We have shown (Heyvaerts and Norman 2002a) that in the vicinity of the 
neutral surface a(a) varies as: 

4t = K = ^ (A-10) 
a{a) i]al 
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where rj is a constant having the same dimension as a and z/ is a positive exponent, strictly 
smaller than unity and most often equal to (1/2). This implies that the edge boundary layer 
is not infinitely extended. Indeed, as a approaches a*, X becomes close to unity and, in this 
limit, the relation between radius and a as given by Eqs.(A-9) and (A-10) reduces to: 

rdr~ a - 1/u da (A-ll) 

Since is strictly larger than unity, this implies that, as a diverges near a*, r remains 

bounded. This result is similar to our finding (Heyvaerts and Norman 2002a) that the angular 
thickness of a neutral sheet boundary layer remains limited and in fact small. Assume, 
similarly, that the boundary layer thickness remains small in this case, so that r varies only 
little. The boundary layer thickness is then obtained, using Eqs.(A-8) and (A-10), by turning 
Eq.(A-9) into the following differential equation for r(X): 

dr _ V a* pgWrf-' ^ 4 ( 2-(2-T)X r \ 

dX p*V2K\ ^*rl ) \xl{l-XT)l-l) 1 ' 

Integrating, the thickness A* of the edge boundary layer is found. Let us denote by k(v, T) 
the integral 

- 1 (2 - (2 - T)X r ) dX 



kM = v ; ' ' (A-13) 

J\ X"(l— X L )2 2 V 

which is convergent at X = 1 since v is less than unity. From Eq.(A-12) we get: 

_ &(,/, iy. / 2W^-V 4 (AU) 

The total outer radius r* of the jet, when under the pressure P ex t, is obtained by adding 
A* to the extent of the domain occupied by the flux tube at the inner limit of the edge 
boundary layer. Since this flux is close to being a*, this radius is almost given by Eq.(A-l) 
with a = a*. The total poloidal current is treated as if it were constant up to a*, whereas in 
fact it decreases to zero in the boundary layer. Then we find: 

Qo/^dT 2 ( P V2n(a')da' 



r\ = - — - • - — f " — exp 



r-1 fig ' Wo pioIy/W) - m(a>)/a(a>) 



The edge density p* is related to P ext by the boundary condition 

Pext = Q*pI (A-16) 
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The neutral layer radius is related to P ext and to / by the neutral surface's Bennet relation, 

p I 2 = 2Q*plrl (A-17) 

The axial density p is related to I by Eq.(A-2). Eliminating p*, p and r* in favour of 
J, Eq.(A-15) takes the form of a relation between I and P ext , which is best written by 
introducing the dimensionless variables D and i defined by: 

p* = p a 2 D (A-18) 

1 = iT (A ' 19) 

The relation between / and P ex t then translates into the following relation between i and D: 
2Q, f T \n (QbPm 1 ^ ( V2(lla, r- Cl(a')da' 



p a 2 E 3 /2 \EJ W \SlJ \ E J 
which is of the form: 

A o^T) +Al -ITd- = ^ ( A_21 ) 
The coefficients A and Ai are positive and $(i) is: 



/ v^figa, f a * n(a')da' 



ypoaoE^ 



*W = -P ^ 7= (A-22) 



'0«*y s 4 Ooa(o'). 

When the external pressure becomes very small, the second term on the left-hand side of 
Eq.(A-21) becomes negligible if % remains finite and a non- vanishing solution is then found 
for i = i sup , the value of i that causes the function <3>(i) to diverge. If $(i) does not diverge 
for finite i, the solution for % as D becomes small moves towards % — 0. This solution exists 
also when diverges for finite %. The second term on the left hand side of Eq.(A-21) is 
then again negligible with respect to the first one because of the very rapid divergence of 
as % approaches zero. The solution is then close to that of the equation 

h = A °i^> (A - 23) 

We conclude that in this limit of vanishing external pressure the outer boundary layer has 
negligible effect on the relation between the current and the external pressure. Actually, 
Eq.(A-23) is just the non-dimensional form of Eq.(A-l). As a result 1 = and I = I sup 
remain the only possible solutions in the limit of deconfinement, even when the jet is bordered 
by a neutral magnetic surface. 
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B. Fast Critical Point Position 

For conical asymptotics the function |rVa| is bounded from above and below with 
n(a) < \rVa\ < X(a) as {t/ta) tends to infinity (Heyvaerts and Norman 1989). Thus, for 
conical asymptotics (d(r\ Va\)/ds) approaches 0. We assume that there is only one fast 
point in this analysis. For cylindrical asymptotics, (d\rVa\/ds) also approaches and the 
fast point may be at infinity. For parabolic field lines, representable as by z = K(a)r p{ ~ a \ 
say, r|Va| is given by: 

vV(a) +r 2 /z 2 
r|VG| = \K'/K + p>logr\ (B " 1} 

For non-constant p(a), r|Va| approaches 0, as r approaches infinity. If p(a) is a constant 
larger than unity, r|Va| approaches a constant value and d(r\Va\)/ds approaches zero. If 
p(a) = 1, rjz goes to a constant as s approaches infinity, and again d(r\Va\)/ds approaches 
0. In general, d(r\Va\)/ds becomes very small at large distances, allowing a fast critical 
point to be located in the asymptotic region. 
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a, A 

Fig. 1. — The function (aE/Vt) versus flux a for some classical wind and for a poloidal field 
of supposedly dipolar symmetry. The equatorial plane, at a = A, is the only null surface. 
Both a and (ctE/Q) diverge at a null surface. The current I (Eq.(23)) can only take values 
between zero and I sup , the minimum value of (aE/Q). This function is assumed here to have 
an absolute minimum at some a* away from the polar axis. The value of / at infinitely large 
distances could then be zero or I sup (see however sections (4.2) and (4.3)). The flux contained 
in a cylindrical jet about the axis at such large distances must be 27ra* (section 4.1). For 
some winds, the function aE/Q could be monotonically increasing. In this case / could only 
vanish at infinitely large distances. For relativistic winds, the function (a(E — c 2 )/Vt) should 
be substituted for (aE/Vt) and the proper current K (Eq.(50)) for /, with no other change 
in the above analysis 
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Fig. 2. — A near-maximal current jet. In an intermediate asymptotic regime, the circum- 
polar current, 1^, may still be close to the maximal value I sup . The left panel represents 
(aE/fl) (or (a(E — c 2 )/Q) for a relativistic wind), I sup and 1^ (or K sup and resp.), as 
well as the run with a of the enclosed total current 1(a) (resp. K(a)), taking into account the 
polar (a = 0) and equatorial (a = A) current- carrying boundary layers. The wind structure 
(Eq.(33) is represented in the right panel. Each arrow indicates the wind speed and direction 
for a set of equidistant values of the flux variable a. The dashed line represents the magnetic 
surface a = a*. The regions close to the pole and the equator are boundary layers where 
the circum-polar current and return current flow. These bounday layers are not explicitly 
represented. 
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Fig. 3. — A near-maximal current wind. In an intermediate asymptotic regime, the circum- 
polar current, 1^, may still be close to the maximal value I sup . In the present case, this 
maximum value is supposedly reached at the polar axis, a = 0. The left panel represents 
(aE/Q) (or (a(E — c 2 )/Q) for a relativistic wind), I sup and 

-^oo (or K SU p and Koo resp.), as 
well as the run with a of the enclosed total current 1(a) (resp. K(a)), taking into account 
the polar (a = 0) and equatorial (a = A) current- carrying boundary layers. The wind 
structure (Eq.(33)) is represented in the right panel. Each arrow indicates the wind speed 
and direction for a set of equidistant values of the flux variable a. The hatched region is the 
boundary layer about the polar axis. 



